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Abstract 

We develop an event-by-event ideal hydrodynamical framework where initial state density fluctuations are present 
and where we use a similar flow-analysis method as in the experiments to make a one-to-one V2 comparison with the 
measured data. Our studies also show that the participant plane is quite a good approximation for the event plane. 
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1. Introduction 

Hydrodynamical calculations based on averaged initial conditions typically have problems in reproducing the 
experimentally observed centrality- and pr-dependence of elliptic flow coefficient V2(pt) at RHIC, see e.g. Fig 7.5 in 
Ref. H]. Especially, in the most central collisions V2{pt) is clearly underestimated. In this talk I will show that this 
problem can be solved by using ideal event-by-event hydrodynamics fZ|. This model reproduces the measured viipj) 
for all centi-ality classes from 0-5% to 30-40% up to pr ~ 2 GeV. 



2. Event-by-event hydrodynamics framework 

Initial states are generated here with a Monte Carlo Glauber (MCG) model. First we distribute the nucleons 
into nuclei using the standard Woods-Saxon density distribution. In the transverse plane the two colliding nuclei are 
separated by an impact parameter b which is sampled from a distribution dN/db oc b. Nucleons / and j from different 
nuclei collide if 

(xi - XjY + (yi - yjY < , (1) 

where ctnn is the inelastic nucleon-nucleon cross section. For collisions at -^Saw = 200 GeV, considered here, we 
take c^;^,;^, = 42 mb. For simplicity we do not include any nucleon finite size effects. As shown in Fig.[T]in order to 
define centrality classes we slice the distribution of A^part so that each interval contains a certain percentage of total 
events. Thus the impact parameter may vary freely in each centrality class. 

MCG gives only the positions of wounded nucleons, but we must start our hydro with an energy density profile. We 
have chosen to distribute the energy density around wounded nucleons using a 2D Gaussian as a smearing function. 



, , K ^ ( {x-xif + {y-yif 



(2) 



where cr is a free smearing parameter controlling the width of the Gaussian. The overall normalization constant K and 
the initial time tq - 0.17 fm are taken from the EKRT pQCD+saturation model |3|. 

We solve the standard ideal hydrodynamic equations df,T^'^ = numerically. We have neglected the small net- 
baryon density since we are interested only in particle production at mid-rapidity. For the same reason we can assume 
longitudinal boost-invariance, which reduces the numerical problem to (2-Hl)-dimensions. To be able to solve the 
hydrodynamical equations we need an Equation of State (EoS) to relate pressure and energy density. Our choice is 
the EoS from Ref. ID. 
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Figure 1: Our centrality class definition for Au+Au collisions at 
y%jv = 200 GeV in terms of the number of participants. From 

m 



Figure 2: Positive-pion pj spectrum for Au+Au collisions at 
yJivjv = 200 GeV calculated with averaged and fluctuating ini- 
tial states varying the Gaussian smearing width a. Data are 
from PHENIX |7|. From fT\. 



We use the conventional Cooper-Frye method to calculate thermal spectra of hadrons emitted from a constant- 
temperature surface T^ec - 160 MeV. To get the individual final state particles for the flow analysis, we sample the 
thermal spectrum. We must also take into account all strong and electromagnetic decays before we can compare with 
the data. For this, we let the unstable thermal hadrons decay one by one using PYTHIA 6.4 |5|. 

Elliptic flow of hadrons is calculated here with respect to two different reference planes. The simplest way is to 
use the reaction plane, defined by the impact parameter and beam axis, as a reference plane. Another way to calculate 
V2 is the event plane method |i6J, where one determines an event flow vector for each event 

Qi = ^(PTi cos(20,), Pn sin(20,)), (3) 

where we sum over every particle in the event and where 0, is measured from the x axis, which is here fixed by the 
impact parameter. The event plane angle i/'2 for each event is then defined to be 

arctan(Q„,y/2„,,) 

•A2 = ^ . (4) 

Then the "observed" elliptic flow, V2{obs) = (cos(2(0/ - 1//2))}, can be calculated with respect to the event plane. 
However, the event plane fluctuates around the "true" event plane due to a finite number of particles from which the 
event plane is determined. We apply the two-subevent method lU to estimate the event plane resolution !??2, which 
corrects V2{obs) for the event plane fluctuations. The final event plane elliptic flow is obtained as 

y2{EP) = y2{obs)/5R2. (5) 



3. Results 

We have plotted in Fig. |2] the transverse momentum spectra of positively charged pions from three hydro calcu- 
lations. We can see that we get more particles at high pr with fluctuating initial states and smaller cr, in which case 
larger pressure gradients are present. Most importantly, we reproduce the pr spectra sufficiently well, so that we can 
meaningfully study the efliptic flow. 
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Figure 3: Elliptic flow of charged particles as a function of pr at different centralities in Au+Au collisions at y%,v = 200 GeV. 
Hydrodynamical results shown are for cr = 0.4 fm. Data are from PHENIX | 8 9 | and STAR 1101 . The statistical errors in the 
experimental data are smaller than the symbol size. 



In Fig. [3]we plot the elliptic flow of charged particles as a function of pr at four different centrality classes. We 
compare the event-by-event calculations V2{RP) and V2{EP) and V2{obs) with V2{RP, ave) which is calculated from an 
averaged initial state. We remind, though, that V2{obs) is a non-physical quantity since it depends e.g. on the number 
of particles (rapidity acceptance) used in the calculation. 

We see that V2{RP) and V2{RP, ave) are very similar at all centralities. This means that fluctuations alone do not 
generate more elliptic flow. Especially this means that the V2 deficit persists in most central collisions even with 
fluctuations. However, we see that V2{EP) fits the measured data surprisingly well at all centralities shown here up to 
Pt ~ 2 GeV. Thus we can conclude that the reference plane definition is indeed very important in V2 calculations and 
thus also in viscosity determination from the measured v'2 data. 

In Fig.|4]we have plotted the elliptic flow divided by the initial eccentricity. In the data and in our V2{EP) points, 
V2{EP) is divided by the participant eccentricity, which is the maximal eccentricity of initial matter distribution. We 
have also plotted V2{RP) divided by the eccentricity calculated with respect to the reaction plane. We notice that both 
calculations are in an agreement with the data and each other 

We have also studied how well the event plane and participant plane are correlated. In Fig.|5]we have plotted the 
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Figure 4: Integrated elliptic flow of charged particles in Au+Au 
collisions at = 200 GeV divided by initial eccentricity. 

The data from PHOBOS QI|. From (2|. 



Figure 5: Correlation of the event plane with the participant 
plane, and with the reaction plane at two different centralities. 
The lines are to guide the eye. From |2|. 



distribution of events as a function of the angle difference from event plane to the participant plane and to the reaction 
plane. We see that the event plane is indeed better correlated with the participant plane than with the reaction plane. 
This suggest, as intuitively expected, that the participant plane is quite a good approximation for the event plane. 
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